Modified Grover's algorithm for an expectation value quantum computer 
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The translation of Grover's search algorithm from its standard version, designed for implemen- 
tation on a single quantum system amenable to projective measurements, into one suitable for an 
ensemble of quantum computers, whose outputs are expectation values of observables, is described 
in detail. A filtering scheme, which effectively determines expectation values on a limited portion 
of the quantum state, is presented and used to locate a single item for searches involving more than 
one marked item. A truncated version of Grover's algorithm, requiring fewer steps than the trans- 
lated standard version but locating marked items just as successfully, is proposed. For quantum 
computational devices which only return expectation values, the truncated version is superior to 
its standard counterpart. This indicates that it is possible to modify quantum algorithms so as to 
reduce the required temporal resources by using the ensemble's spatial resources. 

PACS numbers: 03.67Lx 



I. INTRODUCTION 



The standard formulation of quantum algorithms, such 
as the Deutsch-Jozsa [Q, Shor factorization and 
Grover search algorithms |^ , assumes that the algorithm 
will eventually be implemented on a single quantum sys- 
tem amenable to projective measurements. Currently, 
however, the most advanced realizations of quantum al- 
gorithms use room temperature solution state NMR, 
where an entire ensemble of quantum systems must be 
manipulated and which can only yield expectation values 
of certain observables |, |, | 0, |, ||, [o[ |l|, 
|l5| , [r^ . Consequently various steps of the standard 
algorithm formulation have to be modified for implemen- 
tation on such an expectation value quantum computing 
device. The typical approach has been to translate the 
operations as literally as possible so that the expecta- 
tion value quantum computer effectively mimics its sin- 
gle quantum system relative. In this article we consider 
this process for Grover's algorithm and report a superior 
and significantly different version, which effectively uses 
the ensemble's spatial resources to reduce the required 
temporal resources, for implementation on an expecta- 
tion value quantum computer. 

We conclude this introduction with a description of a 
general scheme for translating the standard version of 
a quantum algorithm into a form suitable for realization 
on an expectation value quantum computer. In section || 
we briefly describe the standard formulation of Grover's 
algorithm. The main result of this article is contained 



in section III which first describes a standard expecta- 
tion value version of Grover's algorithm for searching a 
database containing one marked item and then proposes 
a truncated version requiring fewer steps. Section IV ex- 



tends these ideas to situations where there is more than 
one marked item and proposes a scheme of filtered ex- 
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pectation values for extracting the location of a single 
marked item. Finally in section ^ we briefly consider 
initial attempts to apply these techniques to other quan- 
tum algorithms. 

In the standard formulation of any quantum algorithm 
the first step is to prepare a single L qubit quantum sys- 
tem in a pure state This is followed by applying 
of a sequence of unitary time evolution transformations, 
ill,... , Um, performing |V'/) := f/m • ■ • U2U1 
The time evolution transformations collectively form a 
single algorithm unitary operator, C/aig := Um---U2Ui. 
Finally a computational basis projective measurement 
(PM) is performed on an L' qubit subset of the quantum 
system resulting in a; e {0, ... ,2^ — !}■ The correspond- 
ing collection of measurement operators are the projec- 
tors {Px — \x) {x\ I a; = 0, . . . , 2-^ ~ 1} and the output 

X is produced with probability {ipf \ Px \'4'f) = KsjjV'/)!^- 
Typically this is nonzero for more than one value of x, 
some of which may yield incorrect solutions to the prob- 
lem at hand. However, efficient quantum algorithms have 
the property that the probabilities of such failures are ap- 
propriately small; Shor's algorithm offers an example 

In expectation value quantum computation, such as 
standard NMR implementations, there are two impor- 
tant differences. First, the single quantum system is re- 
placed by an ensemble and in the case of solution state, 
room temperature NMR this is in a highly mixed state. 
However, various pseudo-pure state preparation schemes 
H, |l|, ^ result in a state that, as far as typical NMR 
measurements are concerned, resembles a pure state al- 
though it may lack the entanglement properties associ- 
ated with the pure state ]2l|| . Such preparation schemes 
effectively carry out the same initialization required by 
the single, pure state quantum system formulation. Sec- 
ond, the most basic outputs which an expectation value 
quantum computer can provide are expectation values 
of observables over the entire ensemble; this is the case 
with NMR quantum computers. This stands in contrast 
to situations where individual quantum systems are ac- 
cessible, in which case the basic outputs are the results 
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of projective measurements which can then be used to 
derive expectation values. A translation between the 
PM and expectation value (EV) protocols is essential 
for the development of expectation value quantum com- 
putation. This is accomplished by considering PMs on 
an ensemble of identically prepared quantum systems. 
Firstly a PM can be accomplished in a bitwise fashion 
via projective measurement in the basis {|0) , |1)} for each 
qubit. The resulting outcomes yield the expectation val- 
ues of Pq*^-* = |0) (0|, from which the expectation value of 

tri'^-* = 2Pq'^'' — J^*^) is easily computed (here k labels the 
qubit). The comparative ease with which Pauli operators 
are manipulated versus computational basis projection 
operators results in the standard practice of replacing 
the PM by the EVs of the set of single qubit observables 
{cri'^'' I A: = 1, . . . ,n}. Henceforth we assume that it is 
possible to measure the expectation value of ai'^^ for an 
arbitrary qubit. Finally we note that it is generally as- 
sumed that the algorithm unitary operator (or at least its 
constituents) can be translated straightforwardly into a 
sequence of operations suitable for an expectation value 
quantum computer. This is certainly the case for stan- 
dard NMR quantum computers H |), |§ 

The two protocols are equally powerful for reversible 
deterministic classical computation, in which case \ipf) = 
\s) for some s e {0, . . . , 2" - 1}. The standard EV yields 



(1) 



where s — Sn ■ ■ ■ S2S1 is the bitwise representation of 
the algorithm solution. On the other hand, the PM 
extracts s = s„ . . . S2S1 directly and the equivalence is 
clear. The only requirement is the ability to distinguish 



between (^ai'^'^^ = +1 and (^'^i'"''^ = ^1 for each qubit. 
In quantum computation lijjf) is often a superposition 
of computational basis states, each of which contribute 
to {ipf \ a\t'^ \ipf)- In some cases, such as those presented 
in section IV or ref [|l8|, these contributions cancel and 
inferences such as those offered in Eq. ([l]) no longer ap- 
ply. It is not surprising that the single quantum system, 
pure state versions of algorithms may have to be modified 
so as to be amenable to implementation via expectation 
value quantum computers. Indeed, this modification has 
already been provided for Shor's algorithm pSl]. 
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II. GROVER'S SEARCH ALGORITHM 

Grover's quantum algorithm for searching an unstruc- 
tured database provides a quadratic speedup compared 
to the classical sequential search The rel- 

atively small computational resources required (com- 
pared to Shor's algorithm) has motivated several suc- 
cessful small scale experimental implementations, using 
NMR quantum computers, for databases which conceal 
one marked item in one of as many as eight locations 



1^, 0, ^ |T^, 12, 0. The primary focus in these demon- 
strations appears to have been the successful implemen- 
tation of the preparation and algorithm unitary operator 
steps; here the typical measurement operations aimed 
to verify the accuracy of the aforementioned preceeding 
steps. 

As usual the standard version of the algorithm is for- 
mulated in terms suitable for single quantum systems 
amenable to projective measurements. The problem in- 
volves a database which conceals M marked items in N 
possible locations, where it is assumed that M is known 
(note that there exists a quantum algorithm, which is 
more efficient than its classical counterparts, for finding 
M if it is not known in advance ) . The task is to find 
the location of any one of the marked items, given an 
oracle which can answer the question, "Is a marked item 
located at xl" for any possible location denoted x. The 
database will be represented asX — {0, 1, . . . , A^— 1} and 
the locations of the marked items form a subset, S G X. 
Querying the oracle is equivalent to evaluating 



lix^S 

1 if x e 5* 



(2) 



and it is assumed that this can be done at unit cost for 
any a; G X. A classical sequential search evaluates f{x) 
for a succession of distinct x ^ X until it returns 1 . For 
M <C N this requires approximately N/{M + 1) oracle 
queries on average to be certain of locating at least one 
of the marked items. 

The standard formulation of Grover's algorithm |2^ 
represents database locations in terms of the computa- 
tional basis states, {\x) ,x = 0, 1, . . . ,2^ — 1}, of an L 
qubit data register, where L — 1 < log2(A^) < L. Infor- 
mation about the location of the marked items is supplied 
via an oracle unitary operation, defined on the computa- 
tional basis states as 



Uf\x) := (-1)^(-)|^> 



(3) 



and extended linearly to superpositions of these. The 
database can be extended to one with 2^ locations by re- 
quiring that f{x) := for x = N, . . . ,2^ — 1; henceforth 
we assume that N = 2^ . The data register is initially 
prepared in an unbiased superposition of all possible lo- 
cations, i.e. IV'i) := l/v^2^X]x=o^ 1-^)' which is an un- 
entangled state and typically easy to create. Repeated 
applications of the Grover iterate G := DUf follow the 
preparation; the crucial oracle unitary operations are in- 
terspersed with the unitary "inversion about the average" 




N-l 



2{c))\x) 



(4) 



where (c) = '^^^q Cx/N. Analysis of the algorithm is 



simplified by noting |24, E5| that, after any number of 



applications of G, the state of the data register has the 
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form 



'M 



El 



VN~ M 



(5) 



where a, f3 eM. satisfy a^+/3^ — 1. Thus the data register 
state is conveniently represented by a two dimensional 
real unit vector v := Accordingly the probability 

with which a data register PM will locate a marked item 
is a^. After m applications of G, 

a = sin [(2m + 1)61/2] and /3 = cos [(2m + l)6'/2] (6) 

where cos 6* = 1 — 2M/N, and therefore each invoca- 
tion of the Grover iterate effectively rotates v through 
angle 9. The algorithm may be regarded [^, |2^ as 
a procedure for rotating the undesirable initial state 
Vi = {y^M/N, yj\ — M/N) as close as possible to the 
desired final state v/ = (1,0). If N/M > 1 then 
e w 2yjM/N and thus after l^jN/M 7r/4j applications 
of the Grover iterate v « (1,0), whereupon a PM will 
yield a marked item with probability close to 1. Round- 
ing errors occur when there is no integer m such that 
(2m H- 1)6*72 — n/2 but even in the worst of these cases 
the probability of success is at least 1 — M/N [^ . Each 
application of G queries the oracle once. Thus Grover's 
algorithm provides a quadratic speedup compared to the 
classical sequential search. We shall refer to this rendi- 
tion of the algorithm as the standard PM version. 



III. EV VERSIONS OF GROVER'S 
ALGORITHM: ONE MARKED ITEM 

The first attempt to translate Grover's algorithm to a 
form suitable for implementation on expectation value 
quantum computers results in a standard EV version 
which essentially differs from its standard PM counter- 
part only in the measurement stage. The translation is 
considerably simpler whenever there is only one marked 
item than for the remaining cases (i. e. M > 1), and 
the key ideas behind modifications to the standard EV 
version are illustrated more clearly in the former. Sup- 
pose that there is only one marked item, located at s. 
In the standard EV version the ensemble is initialized 
via a preparation scheme which results in a pseudo-pure 

state corresponding to = l/'V^2^X)^=o^ 1^)- The 
unitary algorithm operator of the standard PM version 
is retained; only the measurement stage must be altered. 
In the standard PM version the data register is approx- 
imately in the state \tpf) — \s) after [iT^/N/'il applica- 
tions of G. The algorithm is essentially deterministic and 



EVs of the single qubit observables, tri'^^ for k — 1 
yield a bitwise representation of s as described earlier. 
Rounding errors result in discrepancies with relative size 
at most 1/A^. Thus if there is only one marked item the 
translation is accomplished by merely replacing the pro- 
jective measurements of the standard PM version with 



expectation values of tri'^'' and determining the binary 
representation of the marked item's location by the sign 

of (ai'''''^ for each data register qubit. This is the stan- 
dard EV version of Grover's algorithm. To date all ex- 
perimental NMR demonstrations of Grover's algorithm 
have taken this approach P, |7|, p|, [lO, 12, p^. 

The standard EV version of the algorithm does not use 
the ensemble's spatial resources (i.e. the replication of the 
quantum computing device via the ensemble members) in 
any essential way. However, by considering the effect of 
terminating the algorithm after significantly fewer than 
[ny/N /4J applications of G, we shall arrive at a trun- 
cated EV version which effectively reduces the required 
temporal resources by taking advantage of these spatial 
resources and which still locates the marked item with 
certainty. It may be surprising that it could be sensi- 
ble to stop after significantly fewer than [tt^/N / 4J appli- 
cations of G since this would leave the data register in 
a superposition of computational basis states somewhat 
distant from the ideal, namely = \s). At the outset 
it is unclear whether, for such superpositions, the sign 

of {^o'z'"''^ could provide any useful information regarding 
the location of the marked item. However, as we shall 
demonstrate below, this is indeed a reasonable strategy. 
In fact, whenever measurements only yield expectation 
values the resulting approach will usually be superior to 
that of the standard EV version. 

As with the standard EV version the system must 
be initialized in the pseudo-pure state corresponding to 
Itpi) — l/v^X^^flo^ 1^)- Then suppose that G is re- 
peatedly applied an arbitrary number of times, resulting 
in the state of the form given by Eq (||). At this stage 
data register EVs give (the discussion is presented for an 
arbitrary number of marked items for later generaliza- 
tion) 



However, 



= -E(-i) 



N - M 



(7) 



Y,{^lf>'=Q => 5](-l)^'==-5](-l)^'= (8) 



and thus 



xifS 



xes 



ai'^)^ii^^^]J2i-ir^. (9) 



M N - M 



xeS 



Finally the normalization requirement, + P'^ — 1, 
yields 



x£S 



where the EV attenuation 



A 



a'^N- M 
N ~M 



(10) 



(11) 
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determines the magnitude of the EV. When there is only 
one marked item and G has been apphed m times 



where the EV attenuation is 

sin^ [{2m +1)9/2] N - 1 



N - 1 



(12) 



(13) 



The binary representation of the marked item's location 
can still be extracted by the usual process of inspecting 
the sign of (ai''^^ for each data register qubit provided 

that Am is sufficiently large. The latter increases mono- 
tonically with respect to m from (for to = 0) to ap- 
proximately 1 in the case of the standard EV version (for 
TO = [f V^J )■ Thus whenever the basic outputs are EVs 
it is entirely feasible to modify the algorithm unitary op- 
erator so that it uses fewer than [f -s/iVJ apphcations of 

G without reducing the probability of successfully locat- 
ing a marked item. The only effect of this truncation is 
to reduce the magnitudes of the EVs by a factor of Am 
where to is the number of applications of G. Therefore we 
propose a truncated EV version of the algorithm which 
uses an L qubit data register initialized in the pseudo- 
pure state corresponding to = l/v'S^X^^^o^ 1^) ^^'^ 
which terminates after the minimum number of appli- 
cations of G such that it is still possible to distinguish 



reliably between ^cri'^-'^ > and {^cr),y ^ < for all data 

register qubits. When there is only one marked item the 
minimum number of applications of G, to, is based on 



the ability to distinguish between ( a 

(fc)\ 



: +Am and 

y = —Am for all qubits. The resolution in measur- 
ing ^ai*'^, and hence the minimum tolerable number of 

applications of G, usually depends on features of the ex- 
perimental setup such as the ensemble size and the signal- 
to-noise ratio of the observed signal. However, in many 
instances it should be possible to use significantly fewer 
invocations of the oracle than required by the standard 
EV version without reducing the probability of locating 
the marked item. In such cases the truncated EV version 
is clearly superior to its standard EV predecessor. 

Proposals for expectation value quantum computing 
devices typically assume that each of the required op- 
erations will be applied simultaneously to all members 
of the ensemble; this is true for existing NMR realiza- 
tions. The truncated EV version of Grover's algorithm 
reduces the required number of constituent operations of 
the unitary algorithm operator, and hence the required 
temporal resources, by effectively exploiting the ensem- 
ble's spatial resources. In this sense too the truncated 
EV version represents a greater departure from the stan- 
dard PM version than the standard EV version does; not 
only are the measurement stages but also the algorithm 
unitary operators different. 



The truncated EV version may appear to violate the 
optimality ^ of the standard formulation of 

Grover's algorithm. However, it will only work if reli- 
able expectation values can be produced. This requires 
an ensemble of quantum systems and it can be argued 
that each oracle invocation, via implementation of G, on 
the ensemble effectively amounts to one oracle query per 
ensemble member. The truncated EV version of Grover's 
algorithm is not intended to substitute the standard PM 
version. Rather in situations, such as NMR, where the 
basic measurement quantities are expectation values (in 
contrast to PMs) it is intended as an improvement on the 
standard EV version. 



IV. EV VERSIONS OF GROVER'S 
ALGORITHM: MORE THAN ONE MARKED 
ITEM 

Translating the standard PM version of the algorithm 
to EV versions is significantly more complicated if there 
is more than one marked item. Here, after [a/ N/M tt/A\ 
applications of G, the data register is approximately in 



the state = 1/ V -^'^ X]a;es 1^) ('^^ henceforth ig- 

nore the deviation from this approximation, which gives 
a probability of M/N of incorrectly locating a marked 
item, since M <C N). A PM will yield a location of one 
of the marked items with certainty. However, merely sub- 
stituting EVs for PMs, as was done for the single marked 
item case, will sometimes not provide enough useful in- 
formation. An extreme example is that where there are 
two marked items, whose locations have binary expan- 



se ... si and s' : 



. s'l where si = 1 



for all k. In this case ^cr^ y = for all fc. If it is known 

that AI = 2 then the only information that this set of 
EVs provides is that s'f. — 1 — Sk- While it is possible to 
determine any single bit for some s e S*, it is impossible 
to infer two or more bits of any s Cz S with certainty. 
In general the strategy of replacing PMs with EVs only 

offers the following limited information: if i^'^^z'^^ ^ 

then Sfc = for some s € S, while (cri'^'''^ < implies 

that Sk ~ 1 for some s G S. However, such bitwise EVs 
ignore crucial inter-qubit correlations, which are effec- 
tively exploited in the PM scheme to yield the location 
of some marked item. The EV outputs for qubits j and 
k only give Sj and sj. for some s, s' G 5*. However, as 
the example above illustrates, there does not necessarily 
exist any single s" € S whose relevant bits have these 
values, (i. e. such that s" = sj and s'l — s'f.). 

Thus we first develop an EV scheme, which exploits 
correlations between qubits, to circumvent this obstacle 
and then use this to provide a truncated EV version for 
Grover's algorithm when there is more than one marked 
item. The key idea, which will eventually allow for a 
bitwise determination of a marked item's location using 
EVs, is a filtered expectation value^ i.e. an EV on one 
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qubit which is conditional on the states of one or more 
of the remaining qubits. A related technique has been 
proposed in connection with a logical labeling scheme for 
pseudo-pure state preparation g. However, to the best 
of our knowledge, application of this to situations such 
as encountered here has not been described explicitly in 
the literature and the version offered here is somewhat 
simpler to that of 

To illustrate the filtered EV technique consider an in- 
stance where ^"'^^''^ ^ 0- Then si = 1 for some s e S". 
The intention is to determine S2 for some s d S such 
that Si = 1. The idea now is to repeat the algorithm 

so that, in the course of measuring (^ai^^^, the final 
data register state is effectively filtered or restricted to 
1^/) = 1/VME.6S,.,=i I^)- Then the sign of (ai^^) 
gives S2, and hence S2S11 for some s € S such that 
si = 1. The restriction could be accomplished by re- 
defining the oracle so as to fix xi = 1. However, this 
can be avoided by comparing the outputs of two runs 
of the algorithm, the second of which includes an addi- 
tional unitary operation. First note that after applying 
?7aig, |V'/) = l/VME.e5l^> and 

(-f)-^ E (-'r- + j^ E i-^r (14) 

Now suppose that [/aig is followed by a correlation opera- 
tion, defined on computational basis states and extended 
linearly to superpositions, 

C2{si = l)\xL...xi) ~ (f5) 

where € {0,f}, and addition is modulo 2. The ar- 
gument of C describes the filtering conditions and the 
subscript the qubit whose EV will be determined sub- 
ject to such conditions. Applying this after the unitary 
algorithm operation gives 

xes 



where 

g{xk ...xi) := XkXk-i ■ ■ ■ X2X1 ® f . (21) 

Averaging the EVs obtained in the two cases where t/aig 
and Cfe+i(sfc = 1,... ,si = 1) t^aig are applied to j-^i) 



after which 

xGS,xi—0 xGS,xi — 1 

(17) 

since ai^Vi^Vi^^ = — cri^^ Averaging the EVs obtained 
immediately after applying the algorithm unitary oper- 
ation, Eq. (^^^and immediately after the correlation 
operation, Eq. (|f7|), gives 

{-'^')-^ E (^lr^ (18) 

xGS,xi — 1 

which is identical to an EV performed on the state 
^/ Vm^xgs xi=i l-^)- Combining the two runs of the al- 
gorithm, one with and the other without the correlation 
operation applied after the algorithm unitary operation 
has effectively filtered out the the desired portion of the 
state. 

When ^o'z^'^ ^ the filtering is conditional on si — 0. 
The only modification to the scheme described above is 
that the correlation operation is 

C2isi=0)\xL...Xi) := {t^x'^Y' \XL.--Xi) (f9) 

or equivalcntly (72(51 = 0) = ci^^ C2{si = l)ai^\ The 
average of the EVs from the two runs will be identical to 
that performed on l/\/ME^gg^^^Q jx). 

This first round of filtered EVs will give S2 and si for 
some s G S. This procedure can be iterated with appro- 
priate modifications of the correlation operation, eventu- 
ally yielding the location of a single marked item. The 
general idea is demonstrated by considering the situation 
where the filtering is conditional ons^ = Sfe_i = . . . si — 
1. The relevant correlation operator is defined via 



(20) 

I 

yields 

E (-1)''"^^ (22) 

XGS^Xh—l;... ,2:1— 1 

and this achieves the desired filtration. Instances where 
any of the filtering conditions is of the form Sj = are 



I 

/ n.^T,\9{xk---xi) 

Cfe+i(sfc = I, . . . , Si = I) |XL . . . Xi) := (^4'=+^^ J \XL... Xi) 
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accommodated by using 



^ X • 

(23) 



Note that g defined in Eq. (]21|) can be computed with 
k muhiphcations and a single additions, which implies 
that the filtering procedure requires 0(log2 N) additional 
operations at worst. 

This filtered EV technique allows an expectation value 
quantum computer to use Grover's algorithm to locate at 
least one marked item with certainty. The additional cost 
is that the algorithm must be run repeatedly and the rel- 
evant correlation operations must be included. However, 
the number of runs of the algorithm is not prohibitive 
since at each level of filtration only one additional run is 
required (to compare with the results from the previous 
level). Therefore at most log2 runs of the algorithm 
are needed to locate a single marked item with certainty. 

This EV version of Grover's algorithm is also amenable 
to the treatment which resulted in the truncated EV ver- 
sion of the algorithm with only one marked item. Again 
suppose that [/aig is modified so that it terminates af- 
ter m applications of G. Applying this to the system in 

the state = l/V2^X]^=o^ 1^) ^^'^ using the filtering 
procedure to select a subset S" C is easily shown to 
result in 



{-If" 



where the EV attenuation is 

sin^ [{2m +1)9/2] N-M 



N-M 



(24) 



(25) 



The only eS'ect of truncating t/aig is a reduction in the 
magnitudes of the EVs by a factor of Am and all the 
filtered EV procedures described in this section for lo- 
cating a marked item still apply. Thus the truncated EV 
version of the algorithm applies for searching a database 
containing any number of marked items and whenever 
the basic measurement quantities are expectation values 
it is superior to the standard (with filtering as necessary) 
EV version. 

The minimum number of applications of G must be 
such that it yields sufficiently large EVs whenever the 
filtering scheme leaves only one marked item. Here it 

is essential to be able to distinguish between (^ai'''''^ = 

+Am/M and (cri''^) = -Am/M. If this is true then it 

will be possible to discern bit values for marked items re- 
liably even when the filtered EV approach leaves more 
than one marked item. In the standard EV version 
Am ~ 1 and in order to locate a marked item with cer- 
tainty the maximum tolerable error is Egtand ■= 1/M. 
Now suppose that for a particular realization of the com- 
puting device, expectation values are determined to an 



accuracy of e < estand- Then it will be possible to suc- 
cessfully infer the location of a marked item after m ap- 
plications of G provided that Am/M > e. Equation ( [25| ) 
implies that this is true if 



1 



1 



m > - arcsin e/e^tand + (1 - e/estand)Af/^ - t;- (26) 
u Z 



li M <s: N then « 2^MlN and the minimum number 
of applications of G for the truncated EV version is 



TOtr 



■ arcsin i/e/estand + (1 - e/cstand)-^/^- 

(27) 



The number of applications of G required for the stan- 
dard EV version to locate a marked item with certainty 

is mgtand ~ 7r/4y^V/M and thus 



TOstand- arcsin \/e/ estand + (1 - e/estand)-^/^- 
TT 

(28) 



Clearly for a given number of marked items, e/estand 
is the determining factor in the reduction in the min- 
imum number of applications of G, which increases 
monotonically, with respect to this parameter, from 
2TOstand arcsin ^ M/N /tt (for e/ estand = 0) to m^tand (for 
e/estand = !)• If e < estand and M ^ N then 



TOtr 



^stand 



ye/fstand + A^/A^. 



(29) 



One remaining issue is the relationship between the en- 
semble's size and it's ability to estimate expectation val- 
ues accurately. Determining expectation values requires 
an ensemble of identically prepared quantum systems, to 
each member of which the same experiment (i.e. the same 
sequence of unitary transformations and projective mea- 
surements) is applied. In this context, expectation values 
are identified with sample averages calculated from out- 
comes of projective measurements, one per member of the 
ensemble. In particular if the experiment is performed on 
an ensemble consisting of n members and the number of 
times that the PMs on qubit k return and 1 are hq 

and ni respectively then the expectation value ^ui'^'' ^ is 

equated with the sample average Xk = {uq — n\)/n. Here 
the distribution is binomial, in which case Chebyshev's 
inequality |M offers 



<e} > l-l/4ne2 (30) 



Prob- 



where e > bounds the error in the estimate. In re- 
alizations such as standard room temperature, solution 
state NMR, n k, 10^° and estimates within very small er- 
ror bounds are easily attained with near certainty. Thus 
the ensemble size is of little concern for standard NMR 
approaches although it may become an issue elsewhere. 



7 



V. OTHER ALGORITHMS 

It remains to ask whether such modifications apply to 
other quantum algorithms. It is unlikely that the trunca- 
tion will be applicable elsewhere as it requires an iterated 
operation, which amongst quantum algorithms is unique 
to Grover's algorithm. On the other hand the filtered EV 
scheme could perhaps be more widely applicable. How- 
ever, initial attempts to apply it to Shor's factorization 
algorithm have been unsuccessful. The quantum part 
of Shor's algorithm seeks to find the period, r, of 
a well defined function of the number to be factorized. 
For the present discussion the essential details are as fol- 
lows. After applying the unitary algorithm operator, a 
computational basis measurement is performed on an L, 
suitably large, qubit data register. The result of this is, 
with high probability, an integer y « 2^k/r for any one 
offc = 0,l,...,r — 1. All possibilities for k are approxi- 
mately equally likely (i. e. « l/?")- An appropriately ter- 
minated continued fractions approximation to gives 
k/r and if k and r are relatively prime, which occurs suffi- 
ciently often, then the denominator of the continued frac- 
tions approximation gives r. It is well known that simply 
substituting EVs for the PMs will not work ||l8| . One pro- 
posed solution is to have the quantum computer perform 
the continued fractions expansion reversibly and follow 
this with EVs [Q. Although these additional steps are 
still computationally efficient, it would be preferable to 
avoid them. One attempt to do so would be to use fil- 
tered EVs to isolate the computational basis states in 
the region of y w 2^fc/r for some desirable choice of k. 
However, although there is a reasonable hope that, by us- 



ing a variant of the binary search, this will succeed with 
an sufficiently small number of additional steps this ap- 
proach is still futile since the amplitudes of the resulting 
EVs will scale as ~ 1/r and hence diminish exponentially 
as the size of r increases. This straightforward approach 
appears to be doomed and no alternatives are clearly vis- 
ible. 



VI. CONCLUSIONS 

In conclusion we have described in detail modifications 
which translate Grover's search algorithm from a formu- 
lation in terms of single quantum systems amenable to 
projective measurements into one suitable for ensemble 
quantum computers where only expectation values are 
available. This included a thorough exposition of a fil- 
tered expectation value technique, which has the poten- 
tial to be applicable more generally. In the context of 
Grover's algorithm, we presented a significantly modified 
version of the algorithm which, for realizations involving 
expectation values, is superior to the standard version. 
The existence of this superior version indicates that it 
may be worthwhile to investigate similar modifications to 
other quantum algorithms although our initial attempts 
to do so have failed. 
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